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^ ; Abstract 

Confined excitonic complexes in two dimensions, consisting of A'^e electrons 
^ • and Nfi holes, are studied by means of Bethe-Goldstone equations. Systems 

. with up to twelve pairs, and asymmetric configurations with N^ 7^ N^ are 

Q ■ considered. The weak confinement regime gives indication of weak binding or 

o : 

even unbinding in the triexciton and the four-exciton system, and binding in 

^ ■ the higher complexes. 
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I. INTRODUCTION 

Although the excitonic matter as a research object has aheady a relatively long history, 
dating back to pioneer theoretical ideas in the sixties |l|] , the study of small excitonic com- 
plexes confined in regions of nanometer scale became possible only recently, and is closely 
related to dramatic advances in semiconductor technology. 

Due to the smallness of the exciton lifetime in semiconductors (hundreds of picoseconds, 
sometimes even greater), signatures of excitonic complexes are looked for mainly in optical 
experiments. Exciton and biexciton peaks in photoluminescence have been identified in 
quantum dots [0,H. Recently, a confocal microscopy technique has been applied to resolve 
the luminescence coming from a single self-assembled quantum dot [§]. The power of the 
excitation laser is such that up to six pairs are excited simultaneously in the dot. Clearly 
distinct lines coming from transitions between multiexcitonic complexes are identified. 

The calculations presented in our paper, although unrealistic, are inspired by the exper- 
imental work IQ. We address a different question, i. e. the dependence of the ground-state 
energy on the number of pairs. From this dependence, the lowest optical absorption line 
related to the creation of an electron-hole pair in the background of A^ photo-excited pairs 
may be extracted. 

On the other hand, small excitonic systems have been widely studied theoretically in 
order to describe the optical properties of bulk and low- dimensional semiconductors. To 
the best of our knowledge, one of the most complete studies is the paper p, in which a 
variational approach is used to compute the properties of two- and tridimensional free (not 
confined) excitons with up to five particles. Our work has also a close connection with that 
paper. 

By using the Bethe-Goldstone (BG) equation to compute the ground-state energies, we 
are able to rise the number of particles in the complex from a maximum of five in the 
variational approach [Q to 24. Asymmetric systems with A'e 7^ A^h are studied as well. In 
the weak confinement regime, we can also address the question about the stability of the free 



complexes. And in fact, our results give indications that the free triexciton and the four- 
exciton system are unstable (or weakly bound), whereas the higher complexes are stable. 
These properties have a distinct trace in the optical absorption, and may be confirmed by 
the experiments. 

The present paper is complementary to [Q, in which a BCS variational approach is 
used to study intermediate-size complexes, containing from 12 to 180 particles. For those 
sizes where both approaches overlap, as a rule BCS gives lower energy values for strong 
confinement, whereas EG gives better energies in the weak confinement regime. 

II. THE BETHE-GOLDSTONE EQUATION 

We will study a two-dimensional model of A'^e electrons and N^ holes confined by a 
parabolic potential in the plane of motion, and interacting via pure Coulomb potentials. 
Only one conduction and one valence band, both ideally parabolic, will be considered, and 
the masses of holes and electrons will be supposed to be equal. By choosing the scales of 



distances and energies as ^hjmiOQ and Tiloq, we get the dimensionless Hamiltonian 

(1) 
where the fi are in-plane coordinates for the particles, rij = \fi — fj\, ujq is the frequency of 



-|(f.|)^|(f4)^^ 



^-^ IT-. . ^-^ r. . ^-^ ^-^ r. . 



the confining potential, k is the dielectric constant of the semiconductor, (3 = JEc/{hujQ), 
and Ec = me^ /{n^fi ) is the characteristic Coulomb energy. 

(3 is the only parameter entering the Hamiltonian (|I]). It may be thought of as the 
inverse of the confinement strength. In the /3 = limit, we have a picture of non- interacting 
fermions which energy levels and one-particle quantum states are those of the 2D harmonic 
oscillator. The electron (hole) states are grouped into shells with "magic numbers" A^g, A^/i = 
2, 6, 12, .... As (3 grows, correlations between particles become more and more important and 
the previous picture is modified. 

One way to go beyond the independent-particle picture is the independent-pair approx- 



iination in which an exact treatment of the two-body correlations is made. Its main com- 
ponent is the Bethe-Goldstone (BG) equation, extensively used in nuclear matter and finite 
nucleus calculations |0Q]. 

The BG equation applies only to fermionic systems. It describes the motion of an 
independent pair of fermions in the system. The rest of the particles excerpt an indirect 
influence on the pair motion through the Pauli principle. The equation takes the form: 

(Ti +T2 + Qa-yV) Ipa^ = E^ylpa^, (2) 

where a and 7 label the states of each fermion in the pair. These states are below the Fermi 
level. The Tj are the one-particle terms in the Hamiltonian, V is the two-body interaction 
potential (Coulomb), and Qa-y is a projection operator given by 

Q^^ = |a7) («7l + E 1/"'^') (/^'^'l ' (3) 

where the sum runs over states above the Fermi level. Q projects a given function onto 
states over the Fermi level. (rir2|a;7) is the unsymmetrized product of two non- interacting 
one-particle eigenf unctions (we use the notation given in 0). Ea^ is the pair energy, and 
ipa-y - its wave function. 

Equation (0) is formally similar to a pair scattering equation, except for the presence of 
the projection operator (and the fact that all of the states in the external quadratic potential 
are bound states). The pair wave function ipa-y is looked for in the form 

^a, = |«7) + E ^;'A' l/^'A') , (4) 

and the total energy is computed from 

^ = E^r + E^«7, (5) 

where e^^ = Ea-y — e^"-* — e^^-* . The corrections coming from the BG equations are proven to 
be equivalent to summing up all the ladder diagrams in the linked-cluster expansion for the 
energy 0. 



Multiplying (0) from the left by (/iA| or (a7|, we get 

(^5.°^ + 4°^ - E^,)C;i + E (Ml l^ l/i'A') C;i, = - (MI V \a^) , (6) 

Ea, = ef + ef + (a7| \/ |a7) + E («7| ^ l/i'A') C;,],,. (7) 

At', A' 

Equation (^ may be seen as a linear system of equations for the coefficients C^Jy from 
which we obtain C"^ = C^^^iEa-y)- Then, the transcendental equation (0) is solved for Ea-y- 

One shall notice that the Coulomb interaction does not change neither the angular mo- 
mentum of the pair nor the spin of the particles. The matrix elements are real and have the 
following properties: 

(a7l V l/xA) = (7«| V |A/x) = (/xA| V \a^) . (8) 

Due to these properties, the matrix entering the linear system is symmetric. 

In order to solve the BG equation we have to identify, given an initial pair, all the possible 
final states above the Fermi level preserving angular momentum and spin. Then we solve 
the linear system and the transcendental equation to get the pair energy. One shall take 
into account that there are three kinds of pairs in the system, namely e-e, h-h and e-h. In 
the second term of (^, a < 7 means that for identical particles (e-e and h-h) a state |q;7) 
should be counted only once. For e-h pairs, however, we should take into account the two 
possibilities, i. e. |ae7h) and |a/j7e). 

In our calculations, up to 272 harmonic-oscillator one-particle states (16 shells) are in- 
cluded. Whenever possible, the e-h, spatial inversion and time-reversal symmetries are used 
to reduce the actual number of equations to solve. For example, in the N^ = Nh = 6 
problem, 15 systems of linear equations with roughly 700 unknowns each, and 15 nonlinear 
equations for the pair energies are solved. 

III. SYMMETRIC {Ne = Nh) SYSTEMS 

In this section we present the results obtained for symmetric systems, where N^ = Nh = 
N. The parameter j3 is varied in the interval (0, 2.5). 



The first step in our Bethe-Goldstone calculation is to define filled and empty levels, 
i. e. the Fermi surface. To this end, a Hartree-Fock calculation was implemented. The 
harmonic-oscillator states were filled in accordance with the Hartree-Fock results. We show 
them in spectroscopic (nuclear) notation in Table |I[ 2p_, for example, means the second 
level with /^ = — 1. Se, and Sh refer to the total electron and hole spins respectively, and L 
is the total angular momentum (along the z-axis). 

A. The biexciton 

For the biexciton, our starting configuration is one in which the first harmonic oscillator 
shell is filled for both electron and holes (see Table |). This state has zero total angular 
momentum and total spin. In order to write down the Bethe-Goldstone equations we shall 
identify all possible pairs below the Fermi level. In the e-e and h-h sectors there is only one 
pair, but we have 4 e-h pairs. For a given initial pair, the number of final pair states above 
the Fermi level depends on the number of harmonic oscillator shells, Ngheiii included in the 
calculations. Figure 1 shows the results for the energy as a function of (3 and N shell- The 
solid line is a two-point Fade approximant, which construction is described below. 

A few remarks shall be given at this point. First, notice that the convergence is slow. 
Second, the EG results do not exactly reproduce the perturbative energies at small values of 
(3. It means simply that the characteristic distances are much smaller than the characteristic 
pair dimensions, and thus the independent pair approximation breaks down. Fortunately, in 
this strong-confinement regime the energy is dominated by the one-particle energies, which 
are properly accounted for by the BG approach. Finally, the EG energies are seem to 
overestimate the actual binding energies at large [3. 

The two-point Fade approximant for the ground-state energy is constructed in the fol- 
lowing way [|r2,|l^]. The asymptotic expansions 



E\p^, = ho + h(3 + hl3^ + 0{P^) (9) 



i?|^_oo = «o/?' + a2 + 0(l//32), (10) 

are used, where 60 = 4, 61 = -2.50662, 62 = -2.92 0, and Oq = -2.1928 [|, 02 = 1. oq 
is the ground-state energy of the free biexciton in two dimensions, and 02 conies from the 
center of mass oscillation, contributing with a one to the ground-state energy. 

The two-point approximant is a rational function interpolating between (0) and (fOl): 



The values of the parameters are the following, pa = 0.525372, p4 = —2.06239, gi = 0.83554 
and g2 = 0.940527. 



B. N=3, 4, 5, 6 and 12 

Next, we present results for A^ = 3,4,5,6 and 12, computed with Ngh^ = 16. The 
occupied harmonic-oscillator states are indicated in Table |. Note that, in some cases, 
the Hartree-Fock approximation does not suggest a unique state, thus we performed BG 
calculations with different occupations, and select the one corresponding to the minimal 
energy. For example, in the A^ = 3 system we performed calculations also for a configuration 
very similar to the one given in Table |l|, but with total momentum L=2. 

The results are drawn in Fig. ^ in a "scaled" form, i. e. Egg/N^^'^ vs. P/N^^^. This 
scaling comes from the dependence 60 ~ 4/3 N^^'^, hi ^ —0.96 A^^/^, 62 ~ —1.65 A^ and 
ao ~ — A^ for large A^ 0, but it is nicely satisfied even for the smallest systems. In this 
scaled drawing, the A^ = 4 cluster clearly distinguishes as the less bound one. 

The results at the largest values of /?, i. e. in the weak confinement regime, can be 
taken as indications of stability or instability of the free clusters. The triexciton seems to be 
unbound (like in three dimensions |]13|), and the four-exciton system - evidently unbound. 
However, the larger clusters are likely to be stable. The situation may be analogous to 
nuclei, where there is a small instability island around atomic number five. 



The differences E^^i — En are shown in Fig. 0. We notice that these differences enter 
the expression 

hu = Egap + El + E,^ + En+1 - En, (12) 

for the frequency of hght creating a new electron-hole pair in the background of A^ pairs. 
Egap is the semiconductor gap, and E^'^ are the electron (hole) confinement energies in the 
z-direction. 

As following from Fig. ^, the absorption peak corresponding to (|T2|) will exhibit an 
interesting behaviour as a function of the number of photoexcited pairs, A^. High values 
of A^ can be reached by rising the laser excitation power. When A^ is around four, we 
shall observe a highly blueshifted line, followed by a redshifted one as the power is further 
increased. The conclusions coming from our oversimplified model have only a qualitative 
predictive power. The effect is so pronounced, however, that we expect it may be observed 
in experiments. 

IV. ASYMMETRIC CASES 

We also studied non-symmetric or non- neutral systems, in which N^ ^ Nh- In particular, 
the cases N^ = 4, A^^ = 2 and Ne = 6, A^^ = 2. The results for their ground state energies are 
shown in Fig. ^ In a non-neutral system, the energy is an increasing function of (3 for low 
f3 values, as for electrons [Q, but for weaker confinement the -/3^ biexcitonic contribution 
dominates over the /3^/^ repulsion. A biexciton plus remaining electrons shall be seen at 
very large (3. The N^ = 6, A^^ = 2 system clearly reveals the change in slope. 

V. CONCLUSIONS 

In the present paper, we have shown that the Bethe-Goldstone equations may be used 
as a powerful method to study small confined excitonic clusters. Pure electronic quantum 
dots may be studied as well, and external electric and magnetic fields can be easily included 

8 



in the calculations. Larger unbalanced systems with only one exciton may be studied to 
obtain the optical absorption and photoluminescence of small electronic quantum dots. 

Our numerical results suggest a small instability island for the free clusters around Ne = 
Nh = 4. We showed that this instability causes the appearance of distinct absorption peaks 
as the laser excitation power is raised. 

More realistic calculations, closely related to the experimental results [Q, are in progress. 
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FIGURES 
FIG. 1. Convergence of the biexciton energy as a function of the number of harmonic-oscillator 

shells included in the calculations. 

FIG. 2. Energies of the N-exciton systems as a function of [3. 

FIG. 3. Energy differences -Eat+i — En for the smallest clusters studied. 

FIG. 4. Energies of the A^'e = 4, iV/j = 2, and Ne = Q, Nh = 2 systems. 
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TABLES 



TABLE L Occupied electron and hole states for the BG calculations. 
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